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Abstract 


The paper introduces a technique to deal with the problem 
of an elastic domain containing an arbitrarily oriented internal 
crack. The problem is formulated as a system of integral equa- 
tions for a fictitious layer of body forces imbedded in the plane 
along a closed smooth curve encircling the original domain. The 
problems of a half plane with a crack in the neighborhood of its 
free boundary and of an infinite strip containing a symmetri cal 1 y 
located internal crack with an arbitrary orientation are consid- 
ered as examples. In each case the stress intensity factors are 
computed and are given as functions of the crack angle. 


1 . INTRODUCTION 

The plane elastostatic problem for a strip of finite width 
containing a crack perpendicular to the boundary was considered 
by Isida [1, 2], and more recently, by Sneddon and Srivastav [3]. 
In [2] and [3] the crack is symmetrically located, whereas in [1] 
its location with respect to the boundaries is arbitrary. In [1] 
and [2] the Kol osov-Muskhel i shvi 1 i functions, and in [3] the 
integral transforms are used to solve the problem in which the 
external load is assumed to be uniform pressure on the crack 
surface . 


This work was supported by The National Science Foundation under 
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The plane problem of two bonded elastic half planes contain- 
ing a crack perpendicular to, and terminating at, the interface 
was recently considered in [4]. The problem of a half plane with 
a crack perpendicular to the boundary is given in [4] as a special 
case. In [4] the singular integral equation governing the problem 
was obtained by using the Mellin Transforms. 

In the solutions given in [1-4] symmetry in the geometry of 
the medium as well as in the external loads with respect to the 
plane of the crack is essential for the formulation of the problem. 
Generally, when the crack problem loses its geometric symmetry it 
becomes necessary to use techniques which are more and more numer- 
ically oriented such as the finite element or boundary collocation 
methods. The technique described in this paper is developed as 
an alternative to the boundary collocation method. The paper con- 
siders the two-dimensional elastostatic problem for a domain 
bounded by a smooth curve and containing a straight crack. The 
problem is formulated by imbedding a fictitious layer of body 
forces in the plane containing the crack along a closed curve 
which encircles the original boundary of the domain. By express- 
ing the boundary conditions one obtains a relatively simple system 
of Fredholm integral equations for the unknown body forces. In 
limit when the outer curve is shrunk on the boundary the integral 
equations become singular. The technique is borrowed directly 
from the potential theory and may be traced back to the works 
of Betti, Somigliana, and Lauricella [5]. Its application to 
problems in elasticity has recently been revived and elaborated 
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in [6-10]*. 


2. FORMULATION OF THE PROBLEM 

Let the smooth curve L be the boundary of a plane elastic 
domain 0, containing a crack along y=0, -a<x<a (Figure 1). Let 
the external loads acting on the body be the tractions (o , t ) 
on L, (concentrated) body forces P., Q. acting at z, (j=l,...,m), 

J J J 

and the tractions cr y , and a~, x x “ y acting on the crack sur- 
faces x+iO and x-iO respectively. The problem is the determina- 
tion of the stress state in the body, particularly the evaluation 
of the stress intensity factors at the crack tips, x = +a, y=0. 

Let t = x + i y be a point on L which is defined by the follow- 
ing parametric equations: 

t = x+iy , x = x ( s ) , y = y(s) , t = t(s) (1) 

where s is the arc length on L measured in the positive direction 
from a fixed point (Figure 1). Consider another smooth curve, 

L q , in the plane enclosing the domain D. It is assumed that L 
and L Q do not intersect. Let the parametric equations of L q be 

*0 = V^o ’ x o =x o< s o>’ *o = V s o>’ t o‘ t o is o ) < 2 > 

where s Q is the arc length on L q . 

Consider now the infinite x-y plane (of the same material) 
containing a crack along -a<x<a, y=0, which is subjected to the 

* 

After the work for the present manuscript was completed the 
application of the same technique to the plane elasticity 
problem for a truncated infinite wedge appeared in [11]. 
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same body forces P . , (j=l,...,m), and the same crack surface 

• ■{■ ■{■ “* f <• 

tractions o yS x xy , o y , x xy , as that acting on D, and a layer of 
unknown body forces, X(t Q ), Y(t Q ) distributed along L . The basic 
idea underlying the method used in this paper, which is directly 
borrowed from the potential theory, is that the body forces X and 
Y can be determined in such a way that the normal and tangential 
stresses, a n , x , on the line L in the infinite plane are equal 
to the specified surface tractions on the boundary L of the orig- 
inal domain D. It is obvious that the stress state within the 
portion bounded by the line L of the infinite plane, which is 

subjected to the external loads P., Q., o\ x + , a~ , x " , X, Y 

J J y xy y xy 

is identical to the stress state in D under the original system 
of external loads. Thus, once X and Y are determined, the problem 
is solved. 


To derive the integral equations for the unknown functions 

X and Y, two elementary solutions for an infinite plane with a 

straight crack are needed. The first is the cracked plane under 

a concentrated body force (X,Y) acting at an arbitrary location 

t Q , and the second is the same plane subjected to the (arbitrary) 

+ + - - 

crack surface tractions a , x , c , x . Let <{>■,, fin and <p ? , fi~ 
be, respectively, the Kol osov-Muskhel i shvi 1 i functions for these 
two problems. The solutions are given in [12]; here only the 
results will be stated. 


4>-j ( z ) = - 


V z > • 


KS <W S . x /_x 
— + Z o + <M Z ) 

z-t„ U-tJ 2 0 


fi^z) = 


9 



*o (z) 


. 2 2 n 1 / 2 ^z-t ^ 1 ( z ) “ 

2(z -a ) ' L o 


Kt 0 )] 


kS 


z- 1, 


1- dd) - I(tj] 


IU) - I (t ) 


0 0 


(2-t 0 ) 2 


2 ' t 0 


( 3 . a-c ) 


X + 1Y 

b " 2 "it(1+k) * 


Id) • (z 2 -a 2 )’/ 2 - z , 


J(z) = —j 


(z‘-a 2 ) 1 / 2 


- 1 . 


(4 . a-c ) 


where X, Y are the components of the body force per unit thickness 

acting at t Q , and k = 3-4v for plane strain and k = (3-v)/(l+v) for 

generalized plane stress, v being the Poisson's ratio. 

a n / .. \ _ / .. \ j i a 

/ 

-a " “ ' -a 

9 n / \ \ J .. 1 3 

/ 

■a " - -a 


a r?\ = C + 1 f R(x)p(x)dx J_ a r q ( x ) dx 

*2 [z) R(zT 2TriR(z) l x-z + 2-rri { x-z 


o ( z ) = , c . + 1 , . f R ( x) p ( x) dx 1_ f q(x) dx 

U 2 (z) RlzJ 2uiR(z) l x-z 2-rri i x-z 

( 5 . a ,b) 


R(z) - (z 2 -a 2 ) 1/2 , 

R(x) * 

R + (x) , 


p(x) * \ l[o+ + a~) - ' 

i(T xy + 

T xy )] ’ 


q(x) * \ E(oJ - o~) - ■ 

i(T xy - 

v> ] • 

( 6 . a-d ) 


where the constant C is determined from the following condition 
of si ngl e- val uedness of displacements [13], 

k J <p 9 (z)dz - / Sl(z’)dl = 0 . (7) 

r r 
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In (7) r is a closed contour encircling the crack. Also, the 
stress states at a point t in the plane corresponding to the two 
basic solutions mentioned above are given in terms of the complex 
potentials <f> k and (k = l,2) by 

0 kx + °ky ■ 2 C*k< t ) + , 

0 ky ' 0 kx + 2iT kxy = (t) - <J» k ( t ) + « k (t)] , 

(k=l ,2) . (8 . a , b) 

In particular, for the infinite plane subjected to uniform pres- 
sure on the crack surface we have a y = CT y = "P 0 > x = T x "y = ° 5 
and 

a 2x + a 2y = 2p o Re ^ t 2_^2jl/2 " ^ * 

a 2 p 0 (t-t) 

a 2y " a 2x + 2l T 2xy = " ( t 2_ a 2^3/2 * (9.a,b) 

At a point t on the line L the normal and tangential 
stresses a n (t) and x n (t) are given in terms of the stress com- 
ponents a x (t), a (t), t (t) b y 

2 ( a - i x ) = ( a + o ) - (o -a +2 it ) e 2 1 01 , (10) 

v n n' v x y' v y x xy' ’ v 

where a is the angle between the outward normal and the x-axis, 

and is a known function of the arc length s (Figure 1). In this 

problem since the crack surface tractions are specified, the 

stress state (a£ x , °2y' T 2xy^ or ^°2n’ x 2n^ assumeci t0 be 
known (see equations 5-8 and 10). 

From (3), (4), (8), and (10) the stress components at a 
point t on L for a concentrated body force S at a point t on L Q 
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may be expressed as 

a ln (t ’ t o ) " iT ln (tjt o ) = sk l^ t > t o ) + Sk 2 (t,t 0 ) (11) 

where k-| and are known functions and are given in the Appendix 
(11) constitutes the Green's function for (a 7 - ix, ). If we 

now assume that the body force S is a continuous function of the 
arc length s Q on L Q , and in the original problem the tractions 
a n (t) and x n (t) on the boundary L of the domain D are known along 
with the crack surface tractions, p, q, and the concentrated body 
forces, P., Q., through superposition we obtain the following 

J J 

integral equation: 


j [S(s 0 )k ] (t,t Q ) + S(s 0 )k 2 (t,t 0 )]ds 0 


m 

+ ^ [R j k l (t - z j) + R j k 2 (t ’ z j )] + [a 2n (t > - iT 2n (t)] 

= o„(t) - i x n ( t) , (t£L) (12) 

where t = t(s), t Q = t 0 (s Q ), and 


R . = 


P. + iQ, 


j 2 tt(1 + k) ’ 


(j=l , . . . ,m) 


(13) 


is the concentrated body force at an internal point z.. (12) 

0 

gives two real integral equations of the first kind for the un- 
known functions X ( s Q ) and Y(s Q ). These are Fredholm-type inte- 
gral equations with bounded and continuous kernels and may be 
solved by the standard numerical techniques (see, for example, 


[14]). 


•fa 

If R. is a distributed body force, then the summation in (12) 
must J be replaced by a line or a surface integral, giving again 
a known function of t. 
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3. 


THE CASE OF L 0 -> L 


From the viewpoint of applications a more effective method 

to solve the problem would be the reduction of the integral 

equation (12) to a singular integral equation by shrinking L Q 

to L. As seen from the Appendix, if t and t are points on the 

same line L the kernel of the integral equation (12) will have a 

Cauchy type singularity. In this problem (12) with the Cauchy 

singularity separated as shown in the Appendix would not be the 

correct integral equation and the limiting process L Q ■+■ L requires 

some care. First, let us consider the infinite plane loaded by 

the crack surface tractions p and q (see (6)), the body forces 

P. and Q-, (j = l,...,m), and the distributed line load S acting on 
J J 

t Q which is now assumed to be a point on the boundary L of the 
given domain D. The left hand side of (12) then gives the ex- 
pression for the stress vector (cr - it ) at any point t inside 
the line L (i.e., in D , see Figure 1). This stress vector must 
approach the known tractions on the boundary, a n (t) -ix n (t), as 
the point t goes to L from inside (i.e., from the positive side). 
In this limiting process the Fredholm kernels K-j(t,t Q ) and 
K 2 (t,t 0 ) given in the Appendix will remain unchanged. The 
remaining part of the integral equation has a Cauchy kernel. 

Hence as t goes to L from inside the limiting process has to be 
treated as taking the boundary value of the holomorphic function 
represented by the Cauchy integral. Thus, by using the Plemelj 
Formulas [13] the integral equation (12) may now be expressed as 
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• 2 i a 

( 1 +k ) ire 1 a S ( s ) + / S(s 0 )(^ r ^ + = J -^)ds () 

L o t -t 

+ / [K 1 (t,t 0 )S(s 0 ) + K 2 (t,t 0 )S(s 0 )]ds 0 

+ | [R j k l< t - z j ) + R j k 2< t ’ z j)] + " ir 2n< t > 

- o n (t) - lT n (t) , (t£L), (14) 

where the bounded kernels (t,t Q ) and k.(-t,z.), ( i =1 , 2 ) are 
given in the Appendix. Referring to the definition of S in 
(4. a), (14) gives two real singular integral equations in the 
unknown functions X(s) and Y(s). 

4. THE HALF PLANE AND THE STRIP 

If the domain D is bounded by simple curves such as a 
circle or infinite straight lines, the analysis given in the 
previous section is somewhat simplified. As an example, consider 
the cracked half plane, x'<d, -°°<y'«x>, shown in Figure 2. If 
the arc lengths s and s 0 are measured from the points B and B q , 
respectively, in the positive direction shown in the figure, 
then t=x+iy and = x o + ly o 9i ve equations of the lines 
L and L Q and may be expressed as 

t = e 1 a ( d + is) , t Q = e ia (d Q +is 0 ) , (15. a, b) 

where d, d Q , and a are constant and are shown in Figure 2. Here 

L and are infinite lines on which the real variables s and s rt 
o o 

vary from to +°°. In the examples considered in this paper, 
namely the half plane and the strip of finite width, it is 
assumed that 
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(j — lj... s (Ti)j ( 1 6 . a , b ) 


° n ( t ) - ix n ( t ) = 0 » Rj = 0 , 

that is, the only external load acting on the medium is the crack 
surface traction. Thus, separating the real and the imaginary 
parts, (12) may be expressed as 

CO 

/ [k n (s 0 ,s)X(s Q ) + k 12 (s 0 ,s)Y(s 0 )]ds = - o 2n (s) > 

— 00 

00 

/ [k 21 (s 0 ,s)X(s 0 ) + k 22 (s 0 ,s)Y(s 0 )]ds 0 = r 2n (s) , 

(_oo<s<oo) t ( 1 7 . a ,b) 

where k.. are real bounded functions obtained from (12), (15) 
and the expressions given in the Appendix, and the input func- 
tions a 2n and x 2n are obtained from (5), (6), (8), (10) and (14) 
in terms of the given crack surface tractions. 

In solving the system of integral equations (16), it is 
expedient to make a change in variables such that the range of 
integration is finite. An appropriate change in variables for 
this purpose is the following: 

s = tan^- , s Q = tan^jp , 

(18) 

- 00 < (s,s Q ) < °° , -1 < (<t>,<P 0 ) < 1 • 

In the case of the strip, if there is no symmetry with 
respect to the y* axis (see Figure 2), there are four unknown 
functions and the system of integral equations giving these 
functions are obtained in a similar way. 

In the alternate formulation of the problem d Q =d, 
t = e ia (d + is), t Q = e 1a (d + is Q ), S = (X+iY)/27r(l+ic) , and instead 
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of the system of Fredhol m- type integral equations, from (14) we 
obtain the following singular integral equation: 


i a 


T- - iY < 5 >] + SfsFT / 


i a 


x ( s o> - 1 Y < s o> 

s o’ s 


ds. 


+ 2?TTh?T 1 ( H i( s * s o )[x(s o) + "(*„>] 
+ H 2 (s,s 0 )[X(s 0 ) - iY(s 0 )]}ds 0 


= - a Q [Re(t/(t 2 -a 2 ) 1/2 
+ x 0 [Im(t/(t 2 -a 2 ) 1 /2 


1 - 


a 2 (t-t)e 2la -, 
2(t 2 -a 2 ) 3/2 J 


) + ie 2la {l 


2(t 2 -a 2 ) 1/2 


t -a* 


(t = e ia (d+is), - 0O< S <°o ). (19) 

where 

H j ( s * s o ) = , ( 3 = 1 ,2) . (20) 

In (19) it is assumed that the crack surface tractions, which are 
the only external loads acting on the half plane, are constant 
and are given by (see Figure 2) 

°y + 1 T xy = ■ ( a o + iT o) 5 ^ y=0, " a<x<a )- (21) 

5. THE STRESS INTENSITY FACTORS 


Let 0 ( z ) and ft(z) be the Kol osov-Muskhel i shvi 1 i functions 
for the infinite plane subjected to the body forces X(s Q ), Y(s Q ), 
P - , Q ., and the given crack surface tractions. From (3-6) it is 

J J 

seen that in the close neighborhood of the crack we can write 


<P(z) 


F-j(z) 

7^1 + 


F 2 (z) 


9 
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Using (3), (15), (22), (23) and (24), the stress intensity 
factors at x=a, y=0 for the half plane shown in Figure 2 may be 
obtained as 

k-, (a) - i k 2 ( a ) = k°(a) - ik°(a) 

00 

+ / [h 1 (a,s )X(s ) + h 2 (a,s 0 )Y(s 0 )]ds 0 , (25) 

where the (Green's) functions h-j and h^ are easily obtained from 
(3), (23) and (24), and k° and k^ are the stress intensity fac- 
tors for the infinite plane subjected to the known crack surface 
tractions (and the body forces P-, Q., if any). Stress intensity 

J J 

factors at x=-a, y=0 may be obtained in a similar way. 

•k 

Here the definition of k] and is such that at y=0, x-a= r 
the asymptotic values of the stresses for small r are given by 

V°x“ k l //Jf • T xy“ k 2 //J7 - 
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6 . 


ON THE SOLUTION OF THE INTEGRAL EQUATIONS 


The solution of the system of Fredholm-type integral equa- 
tions given by (17) is, in principle, stra i ghtf orward . However, 
to improve the effectiveness of the numerical solution, it may 
be necessary to make a change in variables such as. suggested in 
(18) and to use a Gaussi an-type integration formula rather than 
one based on dividing the domain into equal subintervals as 
mentioned in [14]. In this type of problem an appropriate inte- 
gration formula would be the following [15]: 

1 n 

J f (x)dx - l w . f ( x . ) , 

-1 i = 1 1 1 

p n (Xi) = 0 > 

w i = 2/(1 -x i ) 2 [P | ! | (x i )] 2 , ( i = 1 , . . . , n ) ( 26 . a-c ) 

where the Legendre polynomials p n (x) are the related orthogonal 
polynomi al s . 

Two singular integral equations obtained from (19) for the 
unknown functions X ( s ) and Y(s) may be considered as a special 
case of the following system: 

00 M . 00 

Af(x) + dt + / K(x , t)f (t)dt = g(x) , 

-» CO — CO 

( — oo<x<oo) (28) 

where 

A = (a.j) , B = ( b . j ) , K(x,t) = (k i j(x,t)) , 
f(x) = ( f i ( x ) ) , g (x) = ( g .j ( x ) ) , ( i , j = 1 , . . . » n ) . (29) 

The matrices A and B are constant and are such that (A+B) and 
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(A-B) are nonsingular, the kernels k.. and the input functions 

* J 

g.j are known bounded functions, and the functions f. are unknown. 
Generally the quantities given in (29) may be complex. The sys- 
tem of singular integral equations given by (28) may formally be 
regularized as follows: Define 


F ( z ) = 


H(z) = 


± f Mdt 

2iri J t-z at » 

_ 00 

'F(z) . U€S + ) 

.(A+B)” 1 (A-B)F(z) , ( z€S" ) , 


(30. a, b) 


F(z) = ( F j ( z ) ) , 
(A-B)' 1 = C = (c..) , 


H ( z ) = ( H j ( z ) ) , 
-1 


(j 1 , • • . , n ) , 


(A+B)"' = D = (d..) , (31. a, b) 


In terms of the matrix of sectionally holomorphic functions, 
H ( z ) (28) may be expressed as 


oo 

H + ( x ) - H" ( x ) = D[g ( x ) - / K(x , t ) f ( t )dt] , (32) 

- 00 

the solution of which is 


H(z) = J [g(s) - / K(s,t)f(t)dt] ^ + P(z) , (33) 

— 00 — oo 

where P(z) is an arbitrary (matrix) polynomial. Now imposing 
the condition that F(z) should vanish as |z|->-°>, it is clear that 
P(z) must be zero. From 

f(x) = F + (x) - F"(x) = H + ( x ) - C(A+B)H"(x) , (34) 

the solution of (28) may then be obtained as 

oo 

f(x) + / M(x,t)f (t)dt = p(x) , (- °°< x< °°) , (35) 

-00 

where 
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P(x) -SJ£g(x) + ££ /f^d. . 

- oo 

M(x,t) = ( D+C ) K( x , t ) + ^ »*_)- ds , (36. a, b) 

— 00 

M(x,t) = (raj k (x,t)), p ( x ) = ( p j ( x ) ) 5 ( j » k = 1 , . . . , n ) . 

(35) is now a system of ordinary Fredholm integral equations of 
the second kind and may be solved numerically by using the trans- 
formation (18) and the integration formula (26). In (28) if 
K(x,t) = 0, then (35) gives the closed form solution of the sys- 
tem as f(x) = p(x). This is the generalization of the Hilbert 
transforms for one unknown function [16]. 


The foregoing technique would be highly recommended provided 
one can evaluate the kernels m j|<(x,t), ( j ,k = 1 , . . . ,n) in closed 
form. On the other hand if these kernels have to be evaluated 
numerically through the singular integrals given by (36. b), the 


technique could be quite laborious. In this case the following 
simpler and more direct approach may be preferable: Noting that 


/ 


d t 
t-x 


0 


(28) may be expressed as 

°° ,« /. oo 

Af ( x ) + A- / - "tIv X ' dt + / K(x,t)f(t)dt = g(x) , 

- 00 -00 

(_ od<x< o°) . ( 37 ) 

Unlike the solution of the singular integral equations defined 
on arcs, the solution of the singular equations defined on infi- 
nite lines and smooth closed contours (e.g., (28) and (14)) are 
usually bounded and continuous functions. Hence for the purpose 
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of numerical analysis, in these equations the singularity of the 
kernel can always be removed and the integral equation can be 
treated as a Fredholm equation. In (37) note that at t=x the 
integrand in the second term becomes the derivative of f(t) 
which is assumed to be bounded. To solve (37) again a transfor- 
mation such as (18) would be very useful. 

It should perhaps be pointed out that in ordinary applica- 
tions one may prefer to work with the Fredholm equations (12) 
rather than the singular equations (14). One reason for this is 
that the advantage of working with an integral equation of the 
second kind (14) instead of that which is of the first kind (12) 
is somewhat eliminated by the necessity of evaluating the Cauchy 
principal value of singular integrals in (14). Another reason 
for preferring (12) over (14) may be the flexibility it offers 
for improving the accuracy of the results without exhaustive 
numerical work. For example, in solving (12) the results can be 
improved over that obtained through the conventional numerical 
techniques for solving Fredholm-type integral equations either 
by selecting the number of intervals, N on L Q in the integrations 
(and hence, the number of unknowns) less than the number, M of 
the "collocation points", t is (j=l,...,M) and using a least 

square technique, or by selecting N>M and solving the resulting 

★ 

system of equations in some optimal sense . 

The "distance" between L and L q may also be considered as another 
parameter which could be varied to improve the convergence of 
the calculated results. For the straight boundaries a < d Q -d < 2a 
(see Figure 2) appears to give the fastest convergence. Final 
numerical results in this paper were obtained by selecting 
d 0 -d = 1.5a. 
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7. 


NUMERICAL RESULTS 


As a first example we consider an elastic half plane con- 
taining an arbitrarily oriented internal crack near its free 
boundary. The results for uniform normal and shear tractions on 
the crack surfaces are given in Figures 3 and 4. For the special 
cases 0=0 and 0 = j these results reduce to that given in [17] 
and [4], respectively, which are partly reproduced in Figures 5 
and 6. The inclined crack results are given for two values of 
d/a, a and d being the half-length and center-to-boundary dis- 
tance of the crack. Two general observations which could be made 
on the basis of these results are that the stress intensity fac- 
tors increase with decreasing distance of the crack tip to the 
free boundary, and generally for the same crack surface tractions, 
the resistance of the medium to (bri ttl e- type ) fracture would be 
higher for crack angles around 0 = ^ ( i . e . , for cracks nearly 
perpendicular to the boundary) than for 0 = 0 (i.e., for cracks 
parallel to the boundary). 

The results for the infinite strip with a symmetri cal ly 
located internal crack are given in Figure 7. The results are 
obtained for d = 3a and for uniform normal and shear tractions on 
the crack surfaces. In limit when 0 and 0= 0 these results 
reduce to (and agree with) that found in [2,3] and [18], respec- 
tively. In a superposition to obtain the results for other load 
combinations it should be noted that in the results given by 
Figures 3-7 the following crack surface tractions have been used: 

Cy ( x , 0 ) = -o Q , x xy (x,0) = - x Q , ( -a<x<a ) . 
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Figure 3. Normal ( k-| ) and the shear (kg) components of the 
stress intensity factors for a crack in a half 
plane under uniform pressure a . 





Figure 5. Stress intensity factors for a crack parallel and 
for that perpendicular to the free boundary; 
external load: crack surface pressure, a Q . 



Figure 6. Stress intensity factors for a crack parallel and 

perpendi cul ar to the free boundary; external load: 

crack surface shear traction, t =-t . 

xy o 



Figure 7. Normal and shear components of the stress intensity 
factors for a crack in an infinite strip; external 
loads: crack surface pressure, o - - o or crack 
surface shear traction t... = -*r_. y 
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APPENDIX 


The expressions for the kernels k-j and k 2 
k l(t,t 0 ) = f 1 (t > t Q ) + f 2 ( t , t Q ) - f 3 (t s t 0 )e 2la , 

k 2 (t,t 0 ) = f ] (t,t o ) + f 2 (t,t 0 ) - f 4 (t,t 0 )e 2la , 

tTt~ + 2 R( t ) » 

o 


- - 


f o ( t , t ) 


jrItt (— ^ — >Cf (t,t ) + 1 - ] 

2 R (t) t . t() 0 R(t()) 

R(t) = (t 2 -a 2 ) 1/2 , 

0 O L — a 


o t -a t - to 


t -t 


I ( t) - I ( t ) 
f ( t > t o ) ~ tT-~t 


I ( t ) = R(t) - t , J(t) = - 1 , 


(A.l) 
(A. 2) 
(A. 3) 

(A. 4) 
(A. 5) 


+ (t-t)J(t)[^ + J(t 0 )(tV^)} , (A. 6) 

o t - 1 o 


(A. 7) 
( A . 8 . a , b) 


V*.V - zhttt 8^)c-(r|r)(J<t) + J <V - 2f (t.t„» 

l 0 ^0 


+ + ) - j(t ))] 

t -a 

+ 2 - ^j[ + f(t,t 0 ) - Kf(t,t Q )} . (A. 9) 
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As L q shrinks on L (Figure 1), the kernels k-j and k 2 con- 
tain Cauchy type singularity which may be separated as follows: 

f l(t,t 0 ) = t _t + ) , 

o 

f 2 (t,t Q ) = F 2 (t>t o ) ■ 


f 3^ t,t: o^ = t Q -t + F 3^ t#t o^ 



(A.IO.a-d) 


(A. 11 .a-d) 


where F-, (j=l»..,4) and K, , (j=l,2) are bounded functions on L 

J vl 

obtained from (A.3-A.11). 
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